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<^ Abstract 



In this paper, the linear free flexural vibration behaviour of functionally 
graded (FG) size-dependent nanoplates are investigated using the finite ele- 
ment method. The field variables are approximated by non-uniform rational 
B-splines. The size- dependent FG nanoplate is investigated by using Erin- 
gen's differential form of nonlocal elasticity theory. The material properties 
are assumed to vary only in the thickness direction and the effective prop- 
erties for FG nanoplate are computed using Mori-Tanaka homogenization 

^ ■ scheme. The accuracy of the present formulation is tested considering the 

problems for which solutions are available. A detailed numerical study is car- 
ried out to examine the effect of material gradient index, the characteristic 
internal length, the plate thickness, the plate aspect ratio and the boundary 

r> conditions on the global response of FG nanoplate. 
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1. Introduction 



Rapid advancement in the application of micro/nano structures in engineer- 
ing fields, mainly MEMS/NEMS devices, due to their superior mechanical 
properties, rendered a sudden momentum in modelling the structures of nano 
and micro length scale. It has been observed that there is a considerable dif- 
ference in the structural behaviour of material at micro- /nano-scale when 
compared to their bulk counterpart. The difficulty in using the classical the- 
ory is that these classical models fails to capture the size effects. The classical 
model over predicts the response of micro-/nano-structures. Also, in the clas- 
sical model, the particles influence one another by contact forces. Another 
way to capture the size effect is to rely on first principle calculations or molec- 
ular dynamic simulations (MD). But even the MD simulation at micron scale 
is computationally demanding. Hence, there was a need for some theory that 
can incorporate the length scale effect. Several modifications of the classi- 
cal elasticity formulation have been proposed in the literature [ll, y, y, |4|, |5]. 
A common feature of all these theories is that they include one or several 
intrinsic length scales and the particles influence one another by long range 
cohesive forces. In this way, the internal length scale can be considered in 
the constitutive equations simply as a material parameter, also referred to 
as 'nonlocal parameter' and 'internal characteristic length' in the literature. 



Among the various nonlocal theories, Eringen's differential gradient elastic- 
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ity [3[ has received considerable attention to study the static and the dynamic 
characteristics of nano beams |6|, LZj and plates [8|, |9|, [lfj • The theory has been 
applied by several authors to study the axial vibrations [9|, |ll|, ll2l, ll3[ and 
transverse vibrations of nanostructures |l4|, [ll|, |l5|, ll6| • Reddy and Pang 
derived governing equations of motion for Euler-Bernoulli beams and Tim- 
oshenko beams using the nonlocal differential relations of Eringen [3|. The 
classical plate theory (CPT) and first order shear deformation theory (FSDT) 
have been reformulated using the nonlocal differential constitutive relations 



Ill 



[Si, lid, ll2|- It is observed that the increasing effect of the nonlocal pa- 
rameter is to decrease the fundamental frequencies of the structure. In all 
the earlier studies, the nonlocal parameter is kept as a variable and effect of 
the nonlocal parameter on the structural response is studied. One approach 
is to estimate the nonlocal parameter by matching the phonon dispersion 
relations computed by these theories with the lattice dynamics dispersion re- 



lation 18j . Ansari et al. , 19[ used MD simulations to derive the appropriate 
value of the nonlocal parameter. Recently, Eltaher et al., 20] employed Erin- 



gen's differential theory to study the size-effect on the structural response of 



functionally graded nanobeams. Chang [21] studied the axial vibration char- 
acteristics of non-uniform and non-homogeneous nanorods. It was observed 
that the nonlocal parameter has profound impact on the structural response. 

The existing numerical approaches are limited to using radial basis func- 
tions 22f and differential quadrature method [d, [20|, [2l| to study the response 
of nanostructures based on nonlocal elasticity theory. The main objective of 
this paper is to employ Eringen's differential gradient elasticity to study the 
structural response of functionally graded (FG) size dependent nanobeams. 
The plate kinematics is described by FSDT and the governing equations of 
motion are solved within a finite element (FE) framework. The field variables 
are approximated by non- uniform rational B-splines (NURBS) [23(. The ef- 
fect of the nonlocal parameter, the material gradient index, the plate aspect 
ratio, thickness of the plate and boundary conditions on the fundamental 
frequencies are numerically studied. 

The paper is orgranized as follows. A brief over of FG material, Eringen's 
nonlocal elasticity and Reissner-Mindlin plate theory is given in the next 
section. Section [3] presents an overview of NURBS basis functions and spatial 
discretization. The numerical results and parametric studies are given in 
Section HI followed by concluding remarks in the last section. 



2. Theoretical development 



2.1. Functionally graded material 

A rectangular plate (length a, width b and thickness h) with functionally 
graded material (FGM), made by mixing two distinct material phases: a 
metal and ceramic is considered with coordinates x, y along the in-plane 
directions and z along the thickness direction (see Figured]). The material 
on the top surface (z = h/2) of the plate is ceramic and is graded to metal 
at the bottom surface of the plate (z = —h/2) by a power law distribution. 
The homogenized material properties are computed using the Mori-Tanaka 
Scheme 
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25]. 




Figure 1: Coordinate system of a rectangular FGM plate. 



Based on the Mori-Tanaka homogenization method, the effective bulk mod- 
ulus K and shear modulus G of the FGM are evaluated as [24J, |25|, |26|, |27J : 
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Here, Vi (i = c, m) is the volume fraction of the phase material. The sub- 
scripts c and m refer to the ceramic and metal phases, respectively. The 
volume fractions of the ceramic and metal phases are related by V c + V m = 1, 
and V c is expressed as: 



V c {z) 



2z + h 
1h 



n > 



(3) 



where n in Equation 03]) is the volume fraction exponent, also referred to as 
the gradient index. Figure [2] shows the variation of the volume fractions of 
ceramic and metal, respectively, in the thickness direction z for the FGM 
plate. The effective Young's modulus E and Poisson's ratio v can be com- 
puted from the following expressions: 
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Figure 2: Through thickness variation of volume fraction 
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The effective mass density p is given by the rule of mixtures as 



P = PcV c + pmV m (5) 

2.2. Overview of nonlocal elasticity 

Eringen [3|,|29], by including the long range cohesive forces, proposed a nonlo- 
cal theory in which the stress at a point depends on the strains of an extended 
region around that point in the body. Thus, the nonlocal stress tensor a at 
a point x is expressed as: 



a(|x' -x|,r)t(x') dx' (6) 

where t(x) is the classical, macroscopic stress tensor at a point x and the 
kernel function a(|x' — x|) is called the nonlocal modulus, also referred to as 
the attenuation or influence function. The nonlocality effects at a reference 
point x produced by the strains at x and x' are included in the constitutive 
law by this function. The kernel weights the effects of the surrounding stress 
states. From the structure of the nonlocal constitutive equation, given by 
Equation (161), it can be seen that the nonlocal modulus has the dimension 

n i — i 

of (length) 3 |3|, |29[. Typically, the kernel is a function of the Euclidean 
distance between the points x and x'. The nonlocal modulus has a very 
important role to play capturing the nonlocal effects. Ofcourse, ultimately, 
the accurate estimation of a(|x' — x|) would dictate the overall reliability 
of a nonlocal model. Eringen [3|, |30[ numerically determined the functional 



form of the kernel. By appropriate choice of the kernel function, Eringen 



showed that the nonlocal constitutive equation given in integral form (see 
Equation (jSJ)) can be represented in an equivalent differential form as: 

(1 - r 2 L 2 V 2 )<x = t (7) 

where r 2 = -p = (^) , e is a material constant and a and L are the internal 
and external characteristic lengths, respectively. 

2.3. Reissner-Mindlin plate 

Consider a rectangular Cartesian coordinate system (x, y, z) with xy— plane 
coinciding with the undeformed middle plane of the plate and the z- coordi- 
nate is taken positive downwards. The displacement field based on first order 
shear deformation plate theory (FSDT), also referred to as the Mindlin plate 
theory is given by: 



u(x,y,z,t) = u (x,y,t)+z9 x {x,y,t) 
v(x,y,z,t) = v (x,y,t) + z6 y (x,y,t) 
w(x,y,z,t) = w (x,y,t) (8) 



(9) 

The midplane strains e p , bending strain e&, shear strain e s in Equation (0) 
are written as 
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where the subscript 'comma' represents the partial derivative with respect to 
the spatial coordinate succeeding it. The equations of motion are given by: 
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where (In, Ii 2 , I22) = / p(^-,z,z 2 ) dz. The nonlocal membrane stress re- 

-h/2 

sultants N and the bending stress resultants M can be related to the mem- 
brane strains, e p and bending strains e& through the following constitutive 
relations 17| : 



N xx 
N = { N yy } (1 -//V 2 ) = A ( £ /) + B,„.e, ) 

M xx 

M = <( M ra J>(l-//V 2 ) = B /)( s ;) + D b/) s,, (12) 

where the matrices A e = Ay, Bf, e = f?jj and Dbf, = D^; (i,j = 1, 2, 6) are the 
extensional, bending-extensional coupling and bending stiffness coefficients 
and are defined as 

rh/2 _ 

{A tj , Bij, D tJ } = / Q tj {1, z, z 2 } dz, i = 1,2,6 (13) 

J-h/2 

Similarly, the transverse shear force Q = {Q X z,Qyz} is related to the trans- 
verse shear strains e s through the following equation 



Qij = E tJ e s (14) 

where E^ = E = J_ h , 2 QviVj dz; (i,j = 4, 5) is the transverse shear stiffness 
coefficient, Vi,Vj is the transverse shear correction factors for non-uniform 
shear strain distribution through the plate thickness. The stiffness coeffi- 
cients Q i; j are defined as 
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The equations of motion in terms of midplane displacements and rotations 
are obtained by substituting the nonlocal constitutive relations (Equation 



(TT2"|) ) into the equations of motion Equation (ITTJ) . For detailed derivation, 
interested readers are referred to the literature [9|, ll7|, [3l| • 

3. Spatial discretization 

In this study, the finite element model has been developed using NURBS 
basis function. 

Non-uniform Rational B-Splines. The key ingredients in the construction of 
NURBS basis functions are: the knot vector (a non decreasing sequence of 
parameter value, £j < £i+i,i = 0, 1, ■ ■ • , m — 1), the control points, Pi, the 
degree of the curve p and the weight associated to a control point, w. The 



i th B-spline basis function of degree p, denoted by N i>p is defined as [23j, |32| : 

iVi ' 0(O " < else 

AUO = 7^-^*^1(0 + / <+p+1 "/ iV i+ i, P -i(0 (16) 

S«+p — Si Si+p+l — Si+1 

A j9 </l degree NURBS curve is defined as follows: 

m 

C(0 = ^ (17) 

E^, P (0^ 

t=0 

where Pj are the control points and u>j are the associated weights. Fig- 
ure [3] shows the third order non-uniform rational B-splines for a knot vector, 
E = {0, 0, 0, 0, 1/3, 1/3, 1/3, 1/2, 2/3, 1, 1, 1, 1}. The NURBS basis 
functions has the following properties: (i) non-negativity, (ii) partition of 
unity, ^iVjj, = 1; (iii) interpolatory at the end points. These properties, 

i 

make them best suitable to represent the displacement field in elasticity. As 
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Figure 3: non-uniform rational B-splines, order of the curve = 3 
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the same function is also used to represent the geometry, the exact represen- 
tation of the geometry is preserved. It should be noted that the continuity 
of the NURBS functions can be custom tailored to suit the needs of the 



problem. Klassen et al, [33| and Fischer et at, |34| employed NURBS basis 
functions as approximating functions for problems in 2D gradient elasticity. 
The displacement field variable is approximated by the NURBS basis func- 
tions. By multiplying the governing equations of motion by appropriate test 
functions and by employing Bubnov-Galerkin method, we get the following 
set of algebraic equations: 

M(j + K<5 = (18) 

where M is the consistent mass matrix, K is the stiffness matrix and 8 is the 
vector of the degree of freedom associated to the displacement field. After 
substituting the characteristic of the time function 8 = —u 2 8, the following 
algebraic equation is obtained: 

(K-oo 2 M)8 = (19) 

where u is the natural frequency. The frequencies are obtained using the 
standard generalized eigenvalue algorithm. 

4. Numerical Example 

In this section, the free vibration characteristics of FG size- dependent nanoplate 
is numerically studied. Figure H] shows the geometry and boundary conditions 
of the plate. The effect of gradient index n, internal characteristic length //, 
the plate aspect ratio a/b, thickness of the plate h and boundary conditions 
on the natural frequency are numerically studied. The assumed values for 
the parameters are listed in Table [TJ 

In all cases, we present the non-dimensionalized free flexural frequency de- 
fined as: 
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a 



a 



(a) 



(6) 



Figure 4: Plate with boundary conditions: (a) all edges simply supported (SSSS) and (b) 
all edges clamped support (CCCC). 



Table 1: Details of parameters used in the numerical study 



Parameter 



Assumed values 



Gradient index, n 0, 1, 2, 3, 4, 5 

Aspect ratio a/b 1, 2 

Thickness of the plate a/h 10, 20 

Internal length \i (nm 2 ) 0, 1, 2, 3, 4, 5 

Boundary conditions Simply supported and Clamped 
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n = uhj£- (20) 



where u is the natural frequency, p c and G c are the mass density and shear 
modulus of the ceramic, respectively. In the present numerical study, the 
effect of the nonlocal parameter is studied by the frequency ratio, defined as: 

Frequency ratio = — — (21) 

where Qnl and Ql are the non-dimensionalized frequency based on nonlocal 
and local elasticity, respectively. The frequency ratio is a measure of the 
error made by neglecting the nonlocal effects. For example, if -^- = 1, the 
effect of the nonlocal parameter is not significant and becomes important for 
any other value. 

The top surface is ceramic rich and the bottom surface is metal rich. The 
FG plate considered here is made up of silicon nitride (SisN 4 ) and stainless 
steel (SUS304). The mass density p and the Young's modulus E are: p c = 
2370 kg/m 3 , E c = 348.43e 9 N/m 2 for Si 3 N 4 and p m = 8166 kg/m 3 , E m = 
201. 04e 9 N/m 2 for SUS304. Poisson's ratio v is assumed to be constant and 
taken as 0.3 for the current study. Here, the modified shear correction factor 



obtained based on energy equivalence principle as outlined in [35f| is used. 
The boundary conditions for simply supported and clamped cases are (see 
Figure HJ) : 

Simply supported boundary condition: 



u = w = y = on x = 0, a 

v o = w o = 9 x = ony = 0,b (22) 



Clamped boundary condition: 
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u — w — 9 y — v — 9 X — on x = 0, a 

u = w = 9 y = v = 9 X = on y = 0, b (23) 

Before proceeding with the detailed study on the effect of gradient index and 
the characteristic internal length on the natural frequencies, the formulation 
developed is validated against the available results pertaining to the linear 
frequencies for an isotropic nanoplate based on FSDT [36] . Table [2] gives 
the results of non-dimensional frequency for a simply supported square plate 
for different values of nonlocal parameter, the plate thickness and the plate 
aspect ratio. The numerical results from the present formulation are found to 
be in very good agreement with the existing solutions. It can be seen that the 
nonlocal theory predicts smaller values of natural frequency than the local 
elasticity theory. Figure [5] shows the convergence of the first fundamental 
frequency with increasing degrees of freedom (dofs). It can be seen that with 
increasing dofs, the analytical solution is approached and the method shows 
a constant convergence. 

Next, the linear free flexural vibration of FG size-dependent nanoplate is 
numerically studied. The variation of the first natural frequency with the 
gradient index and the nonlocal parameter is graphically shown in Figure [6] 
for a simply supported square plate. 

The values of the nonlocal parameter are assumed to vary between /i = 
(local elasticity) to fi = 5 nm 2 and the gradient index is assumed to vary 
between n = (ceramic) to n = 5. The combined effect of the nonlocal 
parameter /z and the gradient index n is to reduce the natural frequency. 
The decrease in frequency due to increase in the gradient index is due to 
increase in the metallic volume fraction. 

The influence of plate dimension on the frequency ratio for various nonlocal 
parameter is shown in Figure [7] for square plate with plate thickness a/h = 
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Figure 5: The convergence of first fundamental frequency with total degrees of freedom. 
A 3 rd order NURBS curve is used. 



16 



Table 2: Comparison of fundamental frequency uih-J p/G for simply supported plate. 
Tt[36j|. 'Order of curve = 3, 5 x 5 control points. 

a/b a/h \x Ref.tt NURBS* 
0.0930 0.0929 
10 1 0.0850 0.0849 

5 0.0660 0.0659 

1 

0.0239 0.0239 

20 1 0.0218 0.0219 

5 0.0169 0.0170 

0.0589 0.0590 

10 1 0.0556 0.0556 

5 0.0463 0.0464 

2 

0.0150 0.0151 

20 1 0.0141 0.0141 

5 0.0118 0.0118 
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Figure 6: Effect of characteristic internal length and gradient index on the first frequency 
for a simply supported square plate with a/h = 10, a/b = 1 and a =10 nm. 



10 and gradient index n = 5. In this case, the nonlocal parameter is assumed 
to be in the range fi = (0,3). It can be seen that as the length of the 
plate increases, the frequency ratio tend to increase and approach the local 
elasticity solution for considerably larger plate length. From Figure EJ it can 
be seen that the nonlocal parameter has a greater influence on the mode 2 
frequency ratio when compared to the mode 1. Figure [9] shows the influence of 
the nonlocal parameter on the frequency ratio for different plate dimensions 
with a/h = 10 and gradient index n = 5. It can be seen that as internal length 
decreases, the frequency ratio increases and approaches the local elasticity 
solution, irrespective of the plate dimensions. It can be concluded that the 
effect of the nonlocal parameter is prominent for very thin plates. 




20 25 30 35 
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Figure 7: Effect of plate dimension on the frequency ratio (Qnl/Ql for a simply supported 
square plate for various internal length with a/h — 10 and gradient index n= 5. 

The effect of different boundary conditions, viz., all edges simply supported 
and all edges clamped boundary conditions are graphically shown in Fig- 
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Figure 8: Effect of the nonlocal parameter /i on the frequency ratio (SInl/Ql for a simply 
supported square plate for the first two frequencies with a/h = 10 and gradient index n = 
5. 
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Figure 9: Effect of the nonlocal parameter on the frequency ratio (Qnl/Ql for a simply 
supported square plate for various length of the plate with a/h = 10 and gradient index 

n = 5. 
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ure [TD] and tabulated in Table [3] for square and rectangular plate with plate 
thickness a/h = 10 and gradient index n = 5. It can be seen that both the 
boundary condition and the nonlocal parameter has an influence on the fre- 
quency ratio. The higher the nonlocal parameter, the larger is the influence, 
irrespective of the boundary condition. 
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Figure 10: Effect of boundary conditions on the frequency ratio (Qnl/^l for a simply 
supported square and rectangular plate for various internal length with a/h = 10 and 
gradient index n = 5. 



5. Conclusion 

The linear free flexural vibrations of FG size-dependent nanoplate is numeri- 
cally studied using NURBS basis functions within a finite element framework. 
The formulation is based on FSDT and the material is assumed to be graded 
only in the thickness direction according to the power-law distribution in 
terms of volume fraction of its constituents. Numerical experiments have 
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Table 3: Effect of the plate aspect ratio a/b, the thickness of the plate a/h, the nonlocal 
parameter (i and the boundary condition on the natural frequency of a FG plate with a = 
10 and gradient index n = 5. 



a/b a/h \i 



SSSS 



CCCC 



10 



20 



mode 1 mode 2 mode 3 mode 1 mode 2 mode 3 



0.0441 0.1051 0.1051 

1 0.0403 0.0860 0.0860 

2 0.0374 0.0745 0.0746 
4 0.0330 0.0609 0.0610 



0.0758 0.1442 0.1455 

0.0682 0.1157 0.1166 

0.0624 0.0992 0.0999 

0.0542 0.0801 0.0806 



0.0113 0.0278 0.0279 

1 0.0103 0.0228 0.0228 

2 0.0096 0.0197 0.0198 
4 0.0085 0.0161 0.0162 



0.0207 0.0410 0.0421 

0.0186 0.0326 0.0334 

0.0170 0.0279 0.0285 

0.0147 0.0225 0.0229 



10 



20 



0.1055 0.1615 0.2430 

1 0.0863 0.1208 0.1637 

2 0.0748 0.1006 0.1310 
4 0.0612 0.0793 0.0999 



0.1789 0.2251 0.3022 

0.1426 0.1640 0.1960 

0.1218 0.1350 0.1557 

0.0978 0.1051 0.1181 



0.0279 0.0440 0.0701 

1 0.0229 0.0329 0.0464 

2 0.0198 0.0274 0.0371 
4 0.0162 0.0216 0.0283 



0.0534 0.0682 0.0939 

0.0422 0.0491 0.0601 

0.0358 0.0402 0.0476 

0.0287 0.0313 0.0360 
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been conducted to bring out the effect of the gradient index, the nonlocal 
parameter, the plate aspect ratio, thickness of the plate and boundary con- 
dition on the natural frequency of the FG nanoplate. From the detailed 
numerical study, it can be concluded that, the natural frequency decreases 
with increase in the nonlocal parameter and the gradient index. This obser- 
vation is true, irrespective of thick or thin plate, and square or rectangular 
plate considered in the present study. 
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